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The canonical quantization on a single light front is performed for the Abelian gauge fields with 
the Weyl gauge A+ — A~ — coupled with fermion field currents. The analysis is carried separately 
for 1+1 dimensions and for higher dimensions. The Gauss law, implemented weakly as the condition 
on states, selects physical subspace with the Poincare covariance recovered. The perturbative gauge 
field propagators are found with the ML prescription for their spurious poles. The LF Feynman rules 
are found and their equivalence with the usual equal-time perturbation for the S-matrix elements is 
studied for all orders. 
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I. INTRODUCTION 

In the last decade, there is a resurgence of interest in the formulation of quantum field theory on a null-plane 
surface of x + = const, called after Dirac fjl|, a light-front (LF) Though nowadays the basic motivation for it lies 
in the non-perturbative Hamiltonian field theory A3], there are still some problems with a consistent formulation of 
the perturbative LF calculations. In this paper the perturbative approach to the quantum electrodynamics will be 
discussed with a special attention paid to the Abelian gauge vector fields. 

Usually, the LF quantization is performed with the light-cone (LC) gauge condition A_ = A + = being imposed on 
gauge fields, this removes all non-physical gauge modes from the dynamical system and introduces instantaneous (in 
x + ) interactions of currents in the Hamiltonian Q. All these effectively lead to the perturbative gauge field propagator 
with the Cauchy Principal Value (CPV) prescription for the spurious pole at fc_ = 



Oh 
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where the LC gauge vector n^ c = (n+ = l,n^ c = nf c — 0) chooses the LC gauge n LC ■ A = A- = 0. Unfortu- 
nately this result is inconsistent with the equal-time canonical quantization where the LC gauge allows for the 
propagation of non-physical modes and leads to the gauge field propagator with the causal non-covariant distribution 
[k-] AIL = (fc_ +iesgnfc+)~ 1 , so-called the Mandelstam-Leibbrandt (ML) prescription ||. Because the perturbative 
calculations of the gauge-invariant Wilson loops |tJ uniquely indicates the ML prescription as the only consistent 
regularization, then one finds that the LF canonical quantization is perturbatively inconsistent for the LC gauge. 
Recent LF canonical attempts, where one uses two light-fronts x + — const, and x~ — const. 1 as the quantization 
surfaces [||, |)| were only partially succesful. For the free field model, the ML prescription has appeared in the chrono- 
logical product of gauge fields, but one can doubt if this approach can be consistently implemented for interacting 
fields. We believe that another, more fundamental solution to this problem has to be found. 



One possibility is to choose another null axial gauge condition A^ 



A~ 



0, which for the reasons that we will 



explain later, we call the light-front Weyl (LF-Weyl) gauge. In the equal-time formulation both null axial gauges are 
treated on equal footing, therefore we expect the free gauge field propagator in the form 



DZ*« l (x)=i{Q\TA li (x)M0)\Q) 
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where the LF-Weyl gauge vector 



Weyl 



(n 



Weyl 



1, n 



Weyl 



Weyl 

1 <* 



0) chooses the LF-Weyl gauge 



Weyl 



A = A + = 0. However when one LF front x + = const, is fixed as a quantization surface, then the LF- 
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Weyl and the LC gauges are fundamentally different^], the gauge vector n^ c lies on the quantization surface while 
n Weyi j g p er p en( Jicular to this surface. Recently the LF-Weyl gauge has been successfully implemented in the canoni- 
cal formulation of the quantum electrodynamics in the finite LF volume |ll]]. In this model, the canonical DLCQ |l2] ] 
analysis has been carried separately for different subsystems of gauge vector fields and charged fermion fields. This 
simplifying property of the LF Weyl QED is rather a unique situation the LF approach, contrary to the equal-time 
intuitions. We connect it with the simple structure of electromagnetic currents, where no derivatives of fermion fields 
occur. We expect that this simple structure will not be spoiled by the infinite LF volume and in this paper we decided 
to analyse the case of electrodynamics with fermions. We leave the more involved cases of charged scalar fields and 
the non-Abelian gauge interactions for the future publications. 

Our first aim is to study the structure of the gauge field sector therefore we will start the canonical LF quantization 
for the models where the fermion currents are expressed by external arbitrary currents. Because the LF Weyl gauge 
does not fix the gauge symmetry completely, we expect that there will remain dynamical non-physical gauge modes. 
In the canonical equal-time approach [[[3) , |l4| such modes generate the causal LM spurious poles in the propagators 
and introduce positively non-definite Hilbert space of states. We expect that a similar scenario will also appear in the 



LF quantization: the gauge field propagator will take the form (1.2) and the physical states will be chosen by means 
of the Gauss law implemented as a weak condition on states. 

Our second aim is to formulate the LF perturbative Feynman rules for the full QED. As the basic cross-check of the 
LF formulation we take the formal equivalence of the S-matrix elements calculated according to the LF Fynman rules 
with the corresponding elements found within the equal-time approach. We use the same functional method of proof, 
as in p6| , for all orders of perturbation in the coupling constant e. 

This paper is organized as follows. In Section 2 we analyze the 1+1 dimensional model of the Abelian gauge fields 
coupled linearly to external currents. The vector gauge field propagators are given with a consistent definition for 
naively coinciding spurious and physical poles. The higher dimensional model is analysed in Section 3. The infra-red 
(IR) singularities, encountered in the analysis of the independent gauge modes, are dimensionally regularized and 
the final expression for the gauge field propagator is found to have a finite limit in 3+1 dimensions. In Section 4 we 
reintroduce fermion fields and define the perturbative Dyson theory in the interaction representation. The S-matrix 
elements are proved to be (formally) equaivalent to those calculated by the covariant (equal-time) Feynman rules. At 
the end all these results are discussed and the future investigations are outlined. 

All notations and definitions of Green functions are given in Appendix A. In Appendix B the analytical and massive 
regularizations of IR singularities are presented in some details. Appendix C contains calculations for the Poincare 
generators and their commutator relations. 



II. ABELIAN GAUGE FIELDS IN 1+1 DIMENSIONS 

There are three reasons why we decided to start our analysis with the simple 1+1 dimensional gauge field model. 



First, inspecting the ML-pole in (1.2) we find that it contai ns only the longitudinal momenta k±, which are present 
already in 1+1 dimensions. Second, when we put k± = in ( |l.2| ), two poles: physical and spurious coincide at k + = 0. 
If the LF canonical quantization will correctly reproduce these results, then one can also expect its usefullness also for 
more physical models in higher dimensions. Third, the Wilson loop calculations for the Yang-Mills fields in 1 + (D — 1) 
dimensions show singularity at D = 2 |l7|. Though we restrict ourselves to the Abelian fields, it may be instructive 
to compare gauge field sectors in 1 + 1 and higher dimensions. 

In this paper our analysis of the quantum electrodynamics in the 1+1 dimensions will be limited to the canonical 
quantization of Abelian vector gauge field coupled with the external currents jp. These currents describe couplings 
with the fermion fields J M = —eip^^ip when the fermion dynamics is abandoned. This means that no conservation of 
j 71 is supposed, even in opposite, all its components are treated as the arbitrary functions of space-time. 
On can impose the Weyl gauge condition A + = explicitly on the gauge fields and, as the starting point, one can 
take the simple Lagrangian density 



C^y^Ud+A-f + A^r. (2.1) 



1 For example, when the periodic boundary conditions are imposed on the gauge fields, the LC gauge has to be modified jioj 
while the LF-Weyl gauge remains unchanged. 
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This reduced model has only one Euler-Lagrange equation 



which evidently is a dynamical equation and one canonical momentum II = d+A. 
primary constraints, but instead we notice that the Gauss law 



(2.2) 

We stress that there are no 



G 1+1 = d-d+A- + j+ = d_IL + j+ = 
is lost. The canonical quantization is immediate and one obtains the canonical Hamiltonian density 

n\tl = n£M- - c^ yl = i(n) 2 - A_r 

and the nonvanishing canonical commutator 

[H(x + ,x~), A-(x + ,y~)] = -i 6(x~~ - y~). 
They generate the Hamilton equations of motion 

d+U = j- , d+A- = IL~ , 



(2.3) 



(2.4) 



(2.5) 



(2.6) 

which are evidently equivalent to the former Euler-Lagrange equation ( |2.2| ). The Gauss law G 1+1 = cannot be 
imposed strongly as a condition on II, because this would be inconsistent with the above Hamiltonian structure. 
Therefore our quantum theory describes a larger system than the physical electrodynamics and one needs to impose 
the Gauss law rather as a weak condition on physical states 



(phy S '\G 1+1 (x)\phy S ) 



0. 



(2.7) 



The interaction part of fl2.4| ) shows that there is no instant current-current interaction^ quite contrary to the Hamil- 
tonian in the LC gauge | |18| . Therefore the perturbative gauge field propagator is given just by the x + - chro nological 
product of free gauge field operators D__(a;,y) = i(0\T A^(x)A-(y)\0) . When we put j~ = in Eq. ( |2.6|) , the free 
field equations are immediately solved in terms of the cc + -independent operators n(x~) and a_(x _ ) 



Tl(x) = 7r(x ), 
A^{x) = x + ir(x~) + a_(x~), 

which have to satisfy the commutation relation 

[7r(ar),a_(t/")] = -i5(x~ - y~). 

Introducing the Fourier representation for these free fields 

3 +<k-*- t(fc_) + +e - lk - x ~ a(k-)+ 



«-(*-)=/ ^ 
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we get leads to the commutation relations for creation and annihilation operators 

[a(k-),p^(k'_)] = 2in S(k- - k'_) 
[p(k-),at(kL)] =2in S(k--k'_). 



(2.8a) 
(2.8b) 



(2.9) 

(2.10a) 
(2.10b) 



(2.11a) 
(2.11b) 



Because the above expressions show no singularity at the point fc_ = k + = 0, we find that the propagating massless 
gauge fields are less singular than the massless scalar fields in 1+1 dimensions |20|. Also we notice that these gauge 



2 In the full QED, the interaction Hamiltonian contains also dynamical fermions ip + , ijA and has a more complicated structure, 
however with no instantanous interaction ter m. T he perturbative theory based on the full interaction Hamiltonian can be easily 
infered from the results given in the section [VA by omitting the transverse directions and components. 



3 



fields describe non-physical excitations because /r = a* + operators can create Fock states wi th n egative metric. 
The Gupta-Bleuler method can be easily applied here by imposing the weak Gauss law condition fl2.7|) or equivalently 



p(k-)\phys) = 0. (2.12) 

Thus all physical states have zero norm and are created by the operator p^ acting on the vaccun state |0). Usually, 
in the equal-time models, the zero-norm states are also present in the physical subspace while accompanying physical 
photon states. In the LF picture, the physical photons are described by transverse A± gauge fields and such fields 
are evidently absent in 1+1 dimensions, so here we end up with the physical subspace built solely from the zero-norm 
states. 

In the large Hilbert space we immediately find propagators for the independent free fields 

(0|T+A_(x)I%)|0) = Q(x+ - i/+)<0|a_(a:-)7r(tf-)|0) + 6(y+ - a; + )<0|7r(ir)a_(ar)|Q) 

= E 1 F {x-y) (2.13a) 
(0|T+A_(x)A_(y)|0) = x+(0\T 7r(ar)a_(lT)|0) + y + ( Q \ T M^MjT)|0) 

= E 2 F {x-y), (2.13b) 

where the noncovariant Green functions E F (x) and E F {x) are defined in Appendix A and we write a well defined 
Fourier representation for the gauge field propagator 

/oc t2 h, -i 
-oo ( 27r ) [k+ + «e sgn (k-)Y 

Also we notice a property analogous to that for the LC gauge propagators Q 

D—(x + ,x-)=2d-[ d£D F (Z,x-), (2.15) 
Jo 



where D F (x) is the canonical light-front Feynman propagator for a massless scalar field |20) 

f(x+). (2.16) 



17T 

ln|x _ | + — sgn(a; + )sgn(a; _ ) 



Our conclusion is that in 1 + 1 dimensions only propagation of non-physical modes is allowed by the Weyl gauge 
condition and this gives rise to a regular expression for the gauge field propagator. Also we encounter no at fc_ = - 
quite contrary to the expections for a 1 + 1 dimensional model with massless excitations. 

III. CANONICAL QUANTIZATION IN D DIMENSIONS 

In higher dimensions our model has new features, which were absent in the previous case: the vector fields have 
their transverse components Ai and can can propagate also in the transverse directions xj_ . While the first property 
introduces physical photons, then the second one will lead to the unexpected IR singularities. Therefore anticipating 
future problems, we choose to work in D dimensions, where 2 coordinates are LF longitudinal Xl = (x^) and 
remaining d = D — 2 coordinates are LF transverse x±f]. The canonical LF quantization will be carried in D > 4 
dimensions, and the the limit D — > 4 will be taken for all expressions which, being IR-finite, defint the perturbative 
QED in 3+1 dimensions. Here again one can check that the canonical quantization of the vector gauge fields can be 
done independently from the fermion sector, therefore the fermion contribution can be first decribed by the external 
currents j^. As we did in the previous section, we explicitly implement the Weyl gauge condition A + = in the 
Lagrangian density 

C Weyl = d+Ai^Ai - d t A-) + \ (d+A_) 2 - I (diAj - djAi) 2 + A_j~ + Aj. (3.1) 
Our analysis begins with the Euler-Lagrange equations for remaining vector gauge fields 



3 For notation see Appendix A. 
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8 + {d + A--a j A j ) =j-, (3.2a) 
(29+9_ - Al) A % = di (d+A_ - d,A 3 ) + f, (3.2b) 

where = (dj) 2 is the Laplace operator in d = D — 2 dimensions and again the Gauss law G D = is lost, where 

G D = d-(d+A_ - d l A l ) - A d ± A- + j+. (3.3) 

Now we have the primary constraints (in Dirac's nomenclature) connected with the canonical momenta IP = d~Ai — 
diA_, which are not independent canonical variables^,. However these momenta cancel in the canonical Hamiltonian 
density 

n can = n-d+A- + u l d + A, t - c Weyl = in 2 + ^{diAjf + - aj 1 - A-j- (3.4) 

where II~ = d + A_ = H + diAi, so one can proceed freely without any further reference to the constrained momenta 
IP. The remaining fields: Ai,A^ and II are the independent canonical variables with their equations of motion 

(2d+d_ - A ± )Ai = diU + f, (3.5a) 
d+IL = j-, (3.5b) 
d + A-=H + diAi, (3.5c) 

being equivalent to equations ( |3.2|) and the nonzero equal- x + commutators 

[n(x+ > x),A_(x+,y)] = -iS d+1 (x-y) (3.6a) 
[2d_A i (x+,x),A j (x + ,ff)] =-i6 tJ 5 d+1 (x-y) (3.6b) 



where 5 d+1 (x) = 5(x~)5 d (x±). The structure of dynamical equations (3J5) indicates our canonical variables are not 
yet the independent modes, which are very useful in the further formulation. Therefore we propose to introduce the 
following decomposition of vector fields 

A l = C 4 -d i [Ain 1 *(n + 2d i Q), A_ = C_ - MAlp 1 * (IT + 2d 3 C 3 ) (3.7) 

which has the form of gauge transformation and its inverse from C to A has the same functional form. The Green 
function for the Laplace operator is well defined for d > 2 

r A di-i/ T s [d d k ± e^± 1 IW2-1) 
and the last analytical expressin is regular also for noninteger d < 2. 

There are also other regularizations of the inverse Laplace operator strictly in d = 2 dimensions |22| : the analytical 
regularization - when the square pole is changed to 2 — 26 where 1 > 5 > 

and the massive regularization - when the pole at k\ = is shifted by the mass parameter m 

[AJ-^^-^r^H-/ ^W^ = ~h K ^^). (3.10) 



However our decomposition formula ( |3.7[ ) cannot be naively regularized in d = 2 dimensions, because the analytical 
regularization would bring different dimensions for two terms with Cj, while the massive regularization would lead to 
the ill-defined inverse transformation (from d to Ai). The consistent implementation of these two regularizations is 
presented in Appendix [b], while in this sectin we will focus on the dimensional regularization. 
For the independent modes one easily finds nonvanishing commutators 

[U(x+,x),C^(x+,y)] = -i6 d+1 (x-y), (3.11a) 
[2d^C i (x + ,x),C j (x+,y)] =-iS l3 S d+1 (x~y), (3.11b) 
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and effective equations of motion 

(2d+d_ - Al) d = f - 2d l [A d L }- 1 * (8-j- + djf), 

d + C- = -[A d L ]- 1 *(d a l + d-r) 
d+n = r. 



Also the Hamiltonian density (3.4) can be expressed in terms of these new fields 

n eff = ^(diCjf - niAir 1 * [dj + a-r] - a [r - m^ir 1 * P-r + ^)] - c_j- 



(3.12a) 
(3.12b) 
(3.12c) 



(3.13) 



showing again no instantaneous current-current interaction. Therefore, in the interaction picture, the perturbative 
propagators for independent modes will be given by the a; + -chronological products of free fields. The Fourier repre- 
sentation of free fields is defined for positive values of the fc_ momentum 



C_(2) = 
U(x) = 
Ci(x) = 



d d k ± 


f 00 dfc_ r 


(2ir) d 


Jo 


d d k± 


rdL r 


(2ir) d 


Jo ^ 


d d k ± 


r°° dk^ 


(2n) d 


Jo 2tt 2fc_ 



—ik- 



—ik- 



s a(k) 



,+ik-x + 



a T (ft) 



,+ik-x i 



c c,(fc) 



pT(fe) 



>x 4(fc) 



fc-f- — T 



and the commutators for creation and annihilation operators follow from (3.11) 

a{k),p\k')\ = \aHk),p{k')\ = i{2ir) d+1 S d+1 (k-k>) 
Cl (k),c](k')] = (27r) d+1 2/c_ %5 d+1 (fe-fe') 

where <5 d+1 (fc) = 8 d (k 1 _)8(k-f\. The nonzero propagators for the independent modes can be easily found 

(0|T+ C_(f)n(y)|0) = E F (x L - y L )S d (x ± - y x ), 
(0\T+ dMCjivm = S tJ D d F +2 (x - y), 



(3.14a) 
(3.14b) 
(3.14c) 

(3.15a) 
(3.15b) 



(3.16a) 
(3.16b) 



where the covariant Feynman Green function D F +2 (x) is defined in Appendix A3. Next, one can use the decomposi- 
tions (3.7) when calculationf the propagators for the fields A. For the transverse components Ai only the contribution 
of Co fields is taken 



(0\T+M*)Mv)\<>) = (S ik - 29f^[A 



d l-l , 



{5 jl -2d?d?[Ai]- 1 *) D d +\x-y)5 k , 



= 5 l3 D d F +2 (x-y), 



but for the longitudinal component A_ there are more terms 

(Q\T+A_(x)A 3 (y)\0) = ^[Al]- 1 * {2d_D d + 2 (x - y) + E F (x L - y L )5 d { Xl _ -y x )} 
(0\T+A_(x)A_(y)\0) = 2d x _[A d ± ]- 1 * {2dlD d + 2 (x - y) + E F (x L - y L )S d (x ± - y ± )} 



(3.17a) 

(3.17b) 
(3.17c) 



The expression in the above curly brackets can be expressed either as the integral over x + of the covariant Feynman 
Green function 



Q\T + A.(x)A i (y)\0) = d j 
\T+A_(x)A_(y)\0) = 2d- 



dz+D d F +2 (z+ +y+,x-y) 
d^D d+2 (z + + y + ,x~y) 



(3.18a) 
(3.18b) 



4 The commutators ( 3.15a) show that the operators a* and are trivial generalizations of the respective operators in the 1+1 
dimensions and therefore here again they can create negative norm Fock states and the physical subspace has to be chosen by 
means of the weak Gauss law. 
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or as the Fourier integral with two causal poles 



\T+A.(x)A v (y)\0)=i 



d d+2 k e -iHx-y) (k_ n Weyl + k „n™ eyl ) 



(27r) d + 2 k 2 + ie k + +i e' sgn(fc_) 
Finally the Fourier representation can be written in a concise form for all components of field 



■\T+A M (x)A v (y)\0) = i 



d d+2 k e -ik-(x-y) 

(2ir) d + 2 k 2 + ie 



(k^^ 1 + k v n^ 1 ) 
k + + i e' sgn(fc_) 



(3.19) 



(3.20) 



All expressions for propagators have the limit d — > 2 and so we arrive at the canonical LF propagators in 3 + 1 
dimensions 



\T+A^x)A v (y)\0)=i 



d*k 



-ik-(x-y) 



(2tt) 4 k 2 + ie 



(k^^ 1 + k v n^ 1 ) 
k + + i e' sgn(fc_) 



(3.21) 



The relations ( |3.1g| ), which are evidently valid also in 3 + 1 dimensions, can be shown to be the zeroth order (in powers 
of e) terms of the integrated Schwinger-Dyson relations for the full QED in the presence of interactions^]. 
Our above canonical LF analysis gives the same result for the gauge field propagator as the equal-time quantization. 
In both approaches the non-physical modes are allowed to propagate and they introduce causal ML prescription for 
axial spurious poles. 

Before giong to the interacting model with fermions, let us study the Poincare covariance in the gauge field sector. 
We notice that both the choice of a null-surface x + — const as the quantization surface and by the noncovariant 
LF-Weyl gauge condition A + = violate the Poincare covariance. This is quite similar to the equal-time canonical 
procedure for the temporal gauge condition Aq = ]i~3| ], where the Poincare covariance is recovered in the physical 
subspace selected by the weak Gauss law. In our LF formulation we will check this possibility in 3 + 1 dimensions.^ 
First we check that A transforms covariantly 



[M" u , A x {x)] = -i (x"d v - x y d p ) A x (x) - i (g pX A v - g uX A p ) (x). 
Then we analyse the commutator algebra of canonical Poincare generators 

[M^,P X ] = ig vX P p ~ig pX P v + i J d 3 x5 x ~ [g^ G{x)A v (x) - g y -G{x)A p (x)] , 
[M^, M xp ] = +^ [g pX M vp - g vX M pp - g pp M Xv - g vp M Xp ) 

+ i J d 3 i(//+ - x»g p+ ) [g X -G(x)A»(x) - g p -G{x)A x (x)] sym 

-i f d 3 f(xV + -xV + ) [g»-G{x)A»{x)-g v -G{x)A»{x)] sym 



(3.22) 



(3.23) 



(3.24) 



where G = d-(d + A- + diAi) — A±A- is the Gauss law operator in the absence of charged currents and the sym- 
metrization for noncommuting G and A p operators is imposed. In the physical subspace with states selected by the 
weak Gauss law (phys'\G(x)\phys) — one finds that the anomalous terms vanish 



[phys> [g^G{x)A v {x) - g»- G{x)A p {x)} ; 
and one concludes that the Poincare covariance is recovered here. 



phys ) = 



(3.25) 



5 These problems will be discussed in a future publication. 

6 In Appendix D we give the definitions of Poincare generators and present some nontrivial steps of calculations. 
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IV. INTERACTING THEORY WITH FERMIONS 



If one takes the effective Lagrangian density for the gauge sector in the form 

£%l Be = na+A- + d-Aid+Ai - in 2 - ~(^) 2 - nd t A t + A l3 l + A_r, 



(4.1) 



then one can incorporate fermions by resubstituting j M by the fermion currents —eipj^d^ip and adding the kinetic 
terms for fermions. All this leads to the following Lagrangian density for the interacting theory 



£wJ?i = ^d+A- + d-Aid+Ai - ifl 2 - ^{diAjf - TWiAi + iV2iP + d+iP + 



^Weyl —--p-^ — , v 2" o 1 

+ V2ipl (id- - eA-) tp- + + VAC 



where ip± = A±ip, t/>t = ij^ K± and 



iditf + M/3 + ea i A l 



(4.2a) 



(4.2b) 



As usually in the LF formulation, the fermions ip- an d V*— fermions satisfy nondynamical nondynamical equations 

V2{id- -eA-)4i- = £ , -V2 (id- + eAJ) ipl = ft. (4.3) 

For solving these equation uniquely, one has to impose some boundary conditions for the dependent fermion fields. 
Here we choose the simplest and commonly used possibility - the antisymmetric conditions 



lim ip—(x) = — lim ip-(x), lim ij)_(x) = — lim %p_(x), 



x — > — oo 



x -^ — oo 



and this allows us to solve Eqs. (4.3) as 

1 r°° 

^~^ = ~J2j dy ~ ^ id ~ ~ eA -^~ 1 l x ~>y~> x - i -'\€(y~> x - L )> 

dy~ (id- - eA- Y 1 [y~, x~\ x±]£(y~, x±), 



tL(x) 



V2J- 



(4.4) 

(4.5a) 
(4.5b) 



where the Green function (id- — eA-) 1 [x ,y ,x±] is defined by the integral equations 



(id- — eAJ) 1 ' [x ,y ;x±] — (id-) 1 (x — y ) + I dz (id-) 1 (x — z ) eA-(z ,x±) 

J — oo 

x (id- — eAJ) 1 [z~, y~; x±] 

/>oo 

= (id-) (x~ — y~) + / dz~ (id- — eA-) [x~,z~;x±] 



x eA-(z ,x±)(id-) 1 (z —y). 



(4.6a) 



(4.6b) 



It is useful to introduce the matrix notation for the integral operators, by leaving the explicite dependence on coor- 
dinates, denoting the integration over x~ variables by asterisks and writing the integral operators as fractions; for 
example the above equations are written as 



1 



id- — eA- 



1 

id- 



1 

wj 



-eA- * 



1 



1 



1 



* eA-- 



1 

id- 



id- — eA- id- id- — eA 
These equations can be used for generating the series in the arbitrary powers of e, provided the following integral 

dz^ (id) 1 (x~~ — z~) eA-(z, x±) (id-) 1 (z~—y~)^ 



(4.7) 



/OO 
dp-dk- e- ip -^- y ^A-(k-,x±) CPV 
-OO 



CPV 



P- — k- 



(4.8) 
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is a well-defined expression. This means that one has to consider A- fields with their momenta fc_ 0, otherwise in 



Eq. (4.8) one would encounter ill-defined expression ( CPV 



Bearing this limitation in mind we can write 



the infinite series for (id- — eA. 



where 



n—oo , 1 



71=0 s 



id. 



-eA. 



n—oo / 1 



n=0 s 



eA- 



id- 



(4.9a) 



(4.9b) 



(4.9c) 



Thus we write the fermion contribution to the Hamiltonian as 



l f l 
T f er = -j=J d d x x C f * jg- * W-i [a] * £, 



(4.10) 



where the integrand is a local expression in the x± coordinates, and our notation treats £ as a column matrix and 
as a row matrix. 

Thus we end up with the total Hamiltonian, which depends solely on dynamical field variables 



total 



= d 



^Tf + ^diAjf+UdiAi 



1 f°° 1 
-/ Ax^-*W-i[aK, 
V2 J-oc iO- 



while the nonvanishing equal x + (anti)commutators have the canonical LF form 

[ll(x + ,S),A-(x + ,ff)] =-i5 d+ \x-y), 
[2d-Ai(x+, x),A J (x+,y)] = -i6 l3 5 d+1 (x - y), 



(4.11) 



(4.12a) 
(4.12b) 

(4.12c) 



Thus in the LF Weyl gauge A + = the quantum electrodynamics of massive fermion matter fields contains no instan- 
tanous interactions of currents and have all (anti)commutator relations (at equal-x + ) independent of interactions and 
c-numbered. All these properties, while being quite normal in the equal-time formulation, are not such in the front 
description of interacting theories, when the light-front surface probes the fields at the light-like seperated points of 
space-time. 



A. Perturbation theory 



The perturbation theory is defined most easily in the interaction picture, where all quantum field operators have 
free dynamics and the interaction appears in the evolution of quantum states. Thus we also choose to work in this 
representation, however for clarity of further formulas we will omit th e usu ally used subscript i. Free evolution follows 
from a free Hamiltonian, but we notice that the total Hamiltonian (4.11) depends o n the coupling constant e both 
non-locally in the expression W_i[a] and locally in the fields xi and £ T - compare Eq.( 4.2bj ). Thus we define the free 
Hamiltonian Hq as the limit 



H = lim Htotai — 



d d+1 



Itf + ^diAjf + ndiA 



(4.13) 



where now we have 







id l a t + M/3 



(4.14) 



Thus from the (anti)commutator relations ( 4.12| ), which being independent of interactions, remain unchanged, we find 
the free evolution equations 

(4.15a) 
(4.15b) 
(4.15c) 
(4.15d) 

(4.15e) 



(2d+d- 




= sti , 




d+n 


= o , 




d+A_ 


= n + d l A l 


(2d+d- - 


h M 2 ) 4,+ 


= o, 


(2d+d- - 


h M 2 ) $1 


= 0. 



The perturbative propagators are defined as the chronological products of dynamical free fields 
(0 \T+A^x)A u {y)\0) = ig^D F (x - y) + (n^d u + n v d^) [E X F * D d + 2 ) (x - y) , 
T+i> + (x)i>l(y)\o) =iV2A + d x _A d + 2 (x-y,M 2 ) , 



(4.16a) 
(4.16b) 



where for the gauge fields we have repeated our previous analysis, while the result for dynamical fermions is well 
known Jl6| . These expressions will produce all nontrivial contractions in the LF version of the Dyson perturbation 
theory, while the vertices are to be taken from the interaction Hamiltonian, which is defined as the difference between 
the total and the free Hamiltonians Hint — H to tai — Hint- Because all these Hamiltonians are local in the transverse 
directions we can write the expressions for the density (in xj_) 



Hi 



(4.17) 



This expression looks fairly complicated, but fortunately, for our choice of antisymmetric boundary conditions for 
we may consistently introduce the fields and ^_ 



*- = -^=^- * [-idia? - M0\ ip + , *L = i cW + M/3 

y 2 lO— v 2 L 

and further the chronological products for them 



r+*_(a;)Vt(tf)| °) + (° \ T+ ^+( x )^-(y)\ o) = (-d?Y + iM) j°A d + 2 {x - y,M 2 ) 
(o |T*_(a;)* t _(y)| o) = -6^7+7° A£+ 2 (x - y, M 2 ) - 7 + 7°^- * S(x - y). 



(4.18) 

(4.19a) 
(4.19b) 



Now it is easy to check that for the new fermion fields \& = ip+ + \f r _ and ^ = i/A_ + we have a compact expression 
for the propagator 



(0 |Ttt(a;)*(y)| 0) = - (fl^Y* + iM) A F (x - y, M 2 ) - 7+— * S(x - y), 



while a more tedious algebra for the interaction Hamiltonian (4.17) density gives 
Htnt = e* (7"A_ + fA. t ) * 



e 1 



(4.20) 



(4.21) 



These expressions can be compared to tho se obtained in the LC gauge A- = 0. In both cases the chronological product 
f or fr ee fermion fields is the same (4.20) and has one non-covariant non-causal term. Our interaction Hamiltonian 
( 4.21 ) has the infinite number of non-covariant vertices which are generated by W_i[a], contrary to the LC gauge 
Hamiltonian, where there is only one noncovariant term. In the LC gauge case, these two non-covariant terms: from 
the fermion progator and from the interaction Hamiltonian, cancel each other Jl6[ . A similar cancelation should also 
appear in the LF Weyl gauge. However, due to the infinite number of non-covariant terms in (4.21), it is an open 
question if this cancelation will be complete. 
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B. Structure of the S-matrix in the LF perturbative calculations 



Above we have noted that all non-covariant terms are connected with the fermion fields, so in order to prove 
the equivalence of the LF perturbative theory and the usual equal-time formulation we need to study only fermion 
contractions. Keeping the gauge fields as c-numbers, we study Wick's theorem for the chronological product 



Sf er [V, A] = T+ exp -ie / d d x ± * * V * # , 



where 



e 1 

1 + lfi + -Q- * W_i[a] (7 A_ + i>A, 



V = e( 7 -A_+ 7 %) 
Using the Schwinger functional technique [^6| , we derive theone gets the result 



5 /er [*,*;A] = exp[Trln(l-S' F *V)] : exp 



ie / d d xj_^ * V(l — Sf * V)" 1 * * 



where the normal ordering is taken only for fermion fields and 



7+ 1 



S F (x -y) = -i (0 |T+*(aj)* (y)\ 0) = [yfff^ - M) A F (x - y, M 2 ) - * 6(x - y). 



(4.22a) 



(4.22b) 



(4.23) 



(4.24) 



The algebra of Dirac matrices simplifies the contribution for the non-covariant part of the fermion propagator 



7+ 1 
~2~idl 



7+ 1 
~2idl 



* V = — ~ — — — * W-i[o] (7 A- + j j Aj) 



(4.25a) 



and then one easily checks the following factorization property 
1 - S F * V = [1 - S c F ov (7- A. + 7* A,-)] 



e 1 

1 + o7 + -tj- * W_i[o] (7^- +7 fc -4fc) 
2 za_ 



(4.25b) 



where SS™ is the covariant part of the fermion propagator. It is also useful to analyse the non-covariant factor 



Tr In 



e 1 

I + 77V* W_x[a] ( 7 ~A- + 7 fe -4fe) 
2 lO- 



= -trlTrlnW-iIa 1 ]. 



(4.26) 



Gathering all above results we conclude that almost all non-covariant terms have canceled 



Sf er [^,^;A] = : exp 



ie J d d x ± V ( 7 "A_ + j k A k ) * [l - Sf° * (7- A_ + 7 fc A fc )] 1 * * 



exp 



-trlTrlnW_i[a t ] 



exp [Trln [l - Sf" * (7~A_ + 7*A fe )]] 



(4.27) 



So far we have used only algebraic properties, specially when keeping track of the various combinatorical factors. Now 
we may use different arguments connected with the closed loop integrals. We observe that the term TrlnW_i[aT] 
contains the integrals of type 



/ dk 


"1" 




1 




1 


A 


GPV 


k-pi_ 


GPV 


k-p n _ 



iorallpi 7^ 0. 



(4.28) 



GPV 



Thus when fc_ momenta of all insertions of A- are cut at some nozero value, then the non-covariant term is field 
independent may be omitted. In this way we have shown that in the LF perturbative calculation the S-matrix elements 
are the same as those calculated with covariant (equal-time) rules. 
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V. DISCUSSION AND PERSPECTIVES 



In this paper we have presented the canonical quantization procedure, formulated on a single light-front, applied to 
the LF-Weyl gauge QED. Special attention was paid to the free gauge field propagators, where the ML prescription 
arose quite naturally for their spurious poles. This prescription was attributed to the presence of ghostlike modes with 
negative metric. Further implementing the Gauss law weakly as a condition on states, we have separated the physical 
subspace with positive semi-definite metric. In this physical subspace the Poincare covariance was recovered. In the 
analysis of independent modes IR singularities were encountered and several regularizations were studied - all of them 
have lead to the same gauge field propagators in 3 + 1 dimensions. This singularity, specific to the LF formulation, 
appeared for dimensions 2 < D < 4 but not for D = 2 where the gauge field sector is free from IR problems. 
We have noticed similarity between the LF-Weyl gauge in the LF approach and the temporal gauge Aq = in the 
equal-time formulation, which is very encouraging but we are also aware of the limitations. First all this happens for 
the Abelian gauge theory coupled to fermion currents, second there are no self-interactions with derivative couplings 
for gauge fields. Unfortunately, for other types of interaction the light-front quantization would be much more com- 
plicated. Generally one can have nontrivial q-numbered commutators between the vector gauge fields and the matter 
fields and these commutators would depend on interactions (contain coupling constants) so the proper perturbative 
rules could be quite different from those inferred from the free field model. Such things happen in the Abelian model 
with scalar matter fields, where the LF-Weyl gauge leads to the equal- x + commutators between scalar matter fields 
and vector gauge fields which depend on the coupling constant. Similar situation should also occur for the non- Abelian 
case, where the transverse components A\ can have q-numbered commutators with the longitudinal ones At. Such 
phenomena are usually absent in the equal-time approach and one sees that both formalisms may be quite different 
for physically relevant models. 

Though the above observations form a rather narrow limits for the future LF investigations of the LF-Weyl gauge 
models, some possibilities still remain. When the non-Abelian gauge fields are analyzed in the finite volume approach 
(DLCQ) Jl§ the LF-Weyl gauge can be effectively imposed for the zero mode gauge fields. Also it seems worth to 
check the another IR regularization for the light-front procedure, where compactification is introduced only for trans- 
verse space coordinates x± - this, contrary to the common DLCQ method, would leave the parity symmetry x + <-> x~ 
unbroken. For QED one can formally show equivalence of perturbative covariant (equal-time) Feynman rules and the 
LF rules, provided the appropriate regularization for (fc_) _ ™ poles are introduced. Therefore the renormalizability of 
QED in the LF-Weyl gauge should be checked within the 'old-fashioned' Hamiltonian perturbation theory p3| . 



APPENDIX A: NOTATION 



1. Light- front coordinates in D dimensions 

In D dimensions we define 2 longitudinal coordinates x^ — (x ± = x ^ j and take x + as the parameter of dy- 
namical evolution. We denote transverse components x± = (x 2 , . . . , x D ) by the Latin indices Similarly 

we define components of any 4- vector. The metric has nonvanishing components = l,.<?;j = —Sij and the 

scalar product for any 4-vectors decomposes as A ■ B = A + B~ + A~B + — A % B % . Partial derivatives are defined as 
d± = d/dx^, di — d/dx l . Tensor components are defined analogously e.g. T ±A1 = -L(T° M ± T lM ) and summation 
over repeated indices is understood. 

Also we introduce the vector notation for coordinates x — (x~,xi_), which parameterize the light-front surface 
x + = const, and for momenta associated with them k = (fc_,fcj_). The product of such 3- vectors decomposes as 
r\i ' X — - f\i — X " rZ'i Xl . This vector notation is also used for the components of vector gauge field A — (A-, Ai) 



1. Dirac matrices 

The Dirac matrices 7^ satisfy anticommutation relation 

7V+7V = 2<r, (Al) 
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where their components are defined analogously to coordinates e.g. 7 ± = 7 ^? . Thus 7 ± are nilpotent matrices 

(vf = o. 

For the projection operators A± = l/v / 27°7 ± = 1/27 T 7 ± , we have useful relations 



A±A± = A± , A+ + A_ = 1 , 

A+A_ = A_A+ = 7 ±A T = A ±7 T = , 

7 ± A± = A±7 ± 7°A ± = A T 7° , 7 l A± = A±7 l . 

We also use the standard non-relativistic notation pjJ 7 = (3 , 7°7 l = a 1 . 



(A2) 
(A3) 
(A4) 



3. Green functions 

We define the noncovariant Feynman Green functions E f {xl) and E f {xl) as 

0(x+)e- tk - x ~ -&{-x+y k - x 



e x f {x l ) i i 



n 



1 



2n x — ie sgn(x+) 



E%{x L ) i -x+E F (x L ) 



1 



2n x — ie sgn(x+) 

These functions can be also represented by the 2-dimensional Fourier integrals 



E F (x L ) 



E 2 F (x L ) = 1 



d 2 k L 



-ik L -x L 



(2n) 2 k + + ie sgn (fc_) 



d 2 k 7 



-ikr, -XL 



(2tt) 2 [k++ie sgn (£;_)] 2 



and they satisfy obvious relations: E F (0,x ) = 0, 8 + E f (xl) = —E f {xl). 

In D = d + 2 dimensions we have the covariant massive Feynman Green function 



-ik-x 



{2Tr) d + 2 k 2 -m 2 + it 



(A5) 
(A6) 

(A7) 
(A8) 



(A9) 



and its massless limit we denote by D F +2 (x) — lim m _>o A F +2 (x, m2 )- The combination of Green functions, which 
appears in the expression for the longitudinal components of propagator, can be written either as the finite integral 
over x + coordinate or as the Fourier integral 



[2d x _D d F +2 (x)+E F (x L )S d (x 1 _)} = 



00 ,d 



d d ki 



_ x (2ir) d J 2ir 



dk- 



@(x+) 



+ Q{-x+) e 2k - x - 1 I e ik ' s 



A, / dSDp-'fog) 
'0 

d d+2 k e -ik-x 



1 



(27r) d + 2 k 2 +iek + + ie'sgn(fc_) ' 
Similar expressions appear for the massive case in 3+1 dimensions 



[2d*A%(x) + E F (x L )6 2 (x ± )] = (A ± - m 2 

= (Aj_ - m 2 



nr\ I d£A F (Z,x) 




d A k 



-ik-x 



1 



(27r) 4 k 2 — m 2 + ie k + + ie'sgn(k-) 
The analytically regularized inverse Laplace operator in d = 2 dimensions 



(A10) 
(All) 

(A12) 
(A13) 
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evidently is singular when 8 — > 0, however its partial derivative 



<5 



is already finite for 8 — > 

* M -i(,J_* M -i(. x )-,/0^~_£*. CAM) 

In Section 3 many convolutions of various Green functions appear and they are all denoted by * without reference to 
their variables - hopefully leading to no misunderstandings. Here are some explicit examples of them 



[Al]- 1 *n(f)= J d d yi _[Ai}- l (x ± -y ± )U(x-, 



y±) 



[Air 1 * {2d*_D d + 2 (x -y) + Ep(xL - y L )8 d (x ± - y ± )} = J d d w ± [A d L ]- 1 (x ± - w ± ) 
{2d x _D d + 2 (x L - y L , w x - y ± ) + E$.(x L - y L )8 d (w ± - y ± )} . 



APPENDIX B: NONLOCAL AND MASSIVE REGULARIZATIONS 

Strictly in d = 2 dimensions the integral operator di[A±]~ 1 (x±) exists and the parameterization for Ai(x) 

Ai =Ci-di [A J - 1 * (II + TdjCj ) (Bl) 

is regular, but one must be careful not to integrate the partial derivative di by parts because [Aj_] (&x) is ill-defined. 
For field one can introduce the analytical regularizationj] 

A_ = C- - d-[A 2 -2sV * (n + MjCj) (B2) 

by rescaling A- field 

A-{x) -» [Am]- 1 * A_(x) (B3) 

with 1 > 8 > 0, in the Lagrangian density (3.1), but without changing the source term A—j~ . Further, without going 
into details of the canonical procedure, we will present some important points which are different from the respective 
results in Section 3. The canonical commutator for the pair (A-, II) is changed 

[U(x + ,x),A_(x + ,y)] = -tAxIAa-aj]- 1 ^ - yx)S(x~ - y~). (B4) 

and the independent modes vector fields are given by 

Ci = Ai — di[Ax] _1 * (n + 2d j A j ) (B5) 
C- = A- - 9_[A 2 _ 25 ]- X * (n + 2d 3 A 3 ) . (B6) 

Also one commutator for independent modes is changed 

[Tl(x + ,x),C-(x + ,y)] = -iA ± [A 2 - 2S }-\x ± - y ± )8(x- -y~). (B7) 



7 See Appendix A for the definition of [A2-2S] 
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One may worry that these modified commutators violate causality because they do not vanish for spatially separated 
points (in the transverse xj_ coordinates). In the free field case, one can easily check such observables as components 
of the field strength F^ v . First one expresses them in terms of independent modes 



and then finds 



F+- =-A ± [A 2 _ 2S }- 1 *d J C J 
2d-F +j = {S ij A ± -2d i d j )C j . 



[2d_F +t (x+ , x) ,F + _(x+,y)] = -idf ( A ± ) 2 [A 2 _ 25 ] - 1 (x ± - y ± )5(x~ - y~ ) 



(B8) 
(B9) 



(BIO) 



which means that also for the observables causality is violated when 5 ^ 0. Therefore one has the interpretation of 
this analytical regularization as a mere mathematical trick without deeper physical meaning. Keeping this in mind, 
one can proceed further with the canonical procedure and define the Fourier representation for C_ and II free fields 



C-{x) 
U(x) 



oo 



(2tt) 2 



(2tt) 2 



dk_ 
dk_ 

~2tT 



e- lk - s a(k) + e +lk - s a\k) 



-ik-x 



which leads to the commutation relations for the creation and annihilation operators 

a(k),p\k')] = U{k),p(k')] =i(2ir) 3 (kl) 5 S 3 (k- k'), 

where again a'(k) and p'(k') operators may create ghostlike states with negative metric. 
The perturbativc propagators are given by the chronological products of free fields 

(0|T C-(x)Tl(y)\0) = E\{x L - y L ) A ± [A 2 -26}- 1 (x± - y ± ) 
(0\TC i (x)C j (y)\0)=S ij D 4 F (x-y). 



and for the primary gauge fields one obtains the general result 

d A k e - ik < x ~^ 



(0|T A^x)A„(y)\Q) 



(2tt) 4 k 2 +ie 



2 \S 



9^ - (fc±) 



k + + i e' sgn(/c_) 



(Bll) 
(B12) 

(B13) 



(B14) 
(B15) 



(B16) 



where the scaled momentum fc M = (k + ,(k^_) s k^,ki) shows that for 5 > the dimensions of fields A_ and Ai are 



different. Here again one can easily check the integral properties for propagators 

(0|T A^(x)A J (0)\0) =djA ± f dtfAa-a*]- 1 

Jo 

r x+ 

(0|T A_(x)A_(0)\0) = 25_(A ± ) 2 / ^[Aa-atf]" 1 * [A2-25]- 1 * D% (£, x). 

Jo 



(B17) 
(B18) 



The explicit mass term ^rn^A^A^ breaks gauge invariance in the Lagrangian density which no longer describes a 
bona fide gauge theory. However if this mass is treated only as a regularization parameter which will be finally pushed 
to zero for IR finite objects, then one can take the LF-Weyl gauge and start with the Lagrangian density 



C 



Weyl 



d+Ai(d-Ai - diA. 



\ {d+A-f \ {diAj djAif ^Aj + A_j~ + Aif. 



(B19) 



The canonical momenta and commutators for primary canonical variables are the same as in the massless theory 
and the canonical Hamiltonian density has only one extra term +^Af. Now the canonical fields decompose into 
independent modes by the regular invertible transformations 



A i = C i + a i [m 2 -A ± ]- 1 *H. 

A-=C-- <9-[m 2 - AJ- 1 * [2djCj + A ± [m 2 - AJ^n] . 



(B20) 
(B21) 
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The basic difference between the massless and the massive cases is the equation for the longitudinal vector field C_ 



d+C- = m 2 [m 2 - AJ- 1 * II + [to 2 - AJ" 1 * [d t f - A ± [m 2 - AJ" 1 * d-f 
which, for free fields, leads to the Fourier representation with the explicit x + dependence 



n(f) 



00 d 2 k ± 



_ x (2tt) 2 J q 2tT 



C-(x + ,x) = -x + m 2 [m 2 - AJ" 1 *U(x) 



d 2 k ± 



dk- 



(2tt) 2 ,/n 2vr 



-ik- 



x c-(k) + e +tk - s cl(k) 



and one gets nonzero longitudinal components of the propagator 

< 0\T C-(x)C-(y)\0 >= m 2 E 2 F (x L - y L )[m 2 - Ax]" 1 ^ - yx). 

Also the propagator for A_ fields is modified 

< Q|T A_(x)A_(y)\0 >= m 2 E 2 (x L - y L )[m 2 - AJ" 1 ^ - y ± ) 

+2d_A 1 _[m 2 - AJ- 2 * [2dlA 2 + 2 (x - y) + E$(xi - yL)S 2 (x ± - y ± ] 

= 2d*_[ X+ ~ V deAW&a-jft+m* f di \ -'-^ 
Jo Jo Jo 

Finally one derives the Fourier representation of gauge field propagator 



</;" [ drjAp- 2 ^/. .?- ,7). 



<0\TA lt (x)A v {0)\0> =1 



d 4 k 



Ak-x 



(2tt) 4 k 2 - to 2 + it 

2 W>fj,Tlu 



-9^ 



k v n^ L + k^riy 
k + + ie'sgn(k-) 



(B22) 



(B23) 



(B24) 



(B25) 



(B26) 



(B27) 



[k+ + ie'sgn(&:_)] 2 

where the n^riv term is characteristic for the massive regularization and it smoothly vanishes in the limit m 2 — > 



APPENDIX C: POINCARE GENERATORS 

The analysis of the Poincare covariance starts with the definition of the canonical energy-momentum tensor 
ri xv = d u A , d£ + d v A __J?A gl u, c = -p&ff'Ax - g^C 



dd^Ai 
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(CI) 



where we take the Lagrangian density (|3.l|) with the LF-Weyl gauge condition explicitly implemented and we confine 
our discussion to the free Abelian theory in 3+1 dimensions. This canonical tensor may be expressed in terms of the 
symmetric energy-momentum tensor 

= F" X F X V - g^C (C2) 

and one gets 

r itu _ 0M » _ d x (F» x A u ) + 5' 1+ A U G (C3) 
where G is the Gauss law operator in 3+1 dimensions. Similarly we define the canonical angular momentum tensor 



a d^Av d d^Ap 



(C4) 



and then the generators of the Poincare transformations 



1G 



d 2 x ± dx- M +tl ". 



(C5a) 
(C5b) 

(C6) 
(C7) 

where the expected form has regular modifications by the noncovariant terms A^G. These terms arise from the 
noncovariant gauge condition and they generate x + dependence of the Poincare generators M +l 



We notice that the relation (C3) allows to write the expressions for generators as 
P^ = J d 3 x T +p = J d 3 x [6 +AI + A^G] 
AP" = [ d 3 x M +pv = [ d 3 x [(6 +Al + A p, G)x v - (6 +,y + A u G)x"] 



-T—M+ 1 = I d 2 x ± dx- GA i 
dx + 



(C8) 



while other generators are x + independent. For the quantum theory one needs some prescription for noncommuting 
operators and we choose symmetric product for A^G and normal ordering for operators appearing in Q pu . 
After a tedious but straightforward calculation we find all commutators that we need. First for the Poincare generators 
and the independent vector fields A^ = (A_, Aj) | one gets 



P»,A»(x) = -id^A"(x) 
M pu , A x (x) = -i {x"d v - x v d' 1 ) A x (x) - i (g^A" - g vX A^ (x) 



(C9) 
(CIO) 



so A p transforms covariantly. The noncovariant behaviour appears for the field strength tensor F pu , where beside 
the desired relations 



[P»,F xp (x)\ = -id»F xp (x) 
[M» u , F Xp {x)} = -i [x^d" - x u d p ) F Xp (x) 

+ ig vX F pp {x) + ig up F xp (x) - ig pX F vp {x) - ig^F XlJ {x) 
one also gets the anomalous commutator^ 

[M- F-1 {X)] = -i V* - F^x) - iF ij (x) - & ( F+ - {x) + * {X) 



(Cll) 
(C12) 

(C13) 

which is a direct consequence of (|C8f). Fortunately the components F~ J are not present in the Poincare generators 
and one still has covariant commutators 



[M pu , Q +X (x)] = -i (x p d u - x v d») Q +X (x) 

+ ig u+ Q^ x {x) + ig uX Q+^{x) - ig p+ Q vX (x) - ig pX Q +v {x). 

Finally for the (GA P + A P G) term one finds 

M"", (GA P + A p G)(x)} = -i {x p d v - x v d p ) (GA P + A p G){x) 



(C14) 



* [9 



g" + g p -) (GA p + A p G)(x) 



g pp (GA v + A u G){x) - g vp {GA p + A p G){x) 



(C15) 



Gathering all these results one obtains the commutator algebra for Poincare generators (3.23) and (3.24). 



The momentum field II transforms according to the relation II = d+A- — diAi and therefore it will be discarded hereafter. 
9 Here may be any real valued Green function. 
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